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The molecular-crystal model, that describes a one-dimensional electron gas interacting 
with quartic anharmonic lattice vibrations, offers great potentials in the mapping of a rela- 
^^ ■ tively wide range of low-dimensional fermion systems coupled to optical phonons onto quan- 

tum liquids with retarded interactions. Following a non-perturbative approach involving 
non-Gaussian partial functional integrations of lattice degrees of freedom, the exact expres- 
sion of the phonon-mediated two-electron action for this model is derived. With the help 
of Hubbard-Stratonovich transformation the charge-density-wave instability is examined in 
the sequel, with particular emphasis on the effect of the quartic anharmonic phonons on the 
charge-density-wave transition temperature. 
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,,. 1. Introduction 

^^ I Standard electron-phonon models^"^^ have served as a rich workshop for studies of phase 

[^ I transitions in one-dimensional(lD) systems with coexisting structural instabilities and elec- 

^^ . tronic correlations. However, while these standard models provide relevant qualitative insights 

onto physical mechanisms of the phase transitions, potential materials reported to date display 
rather diverse lattice dynamics as well as a great variety of local structures hence bringing 
r> ' about the fundamental question of their universality.^"^-' On another hand, failures of tra- 

^ ■ ditional formalisms such as the Hartree-Fock and weak-coupling expansions of the classic 

electron-phonon interactions, paved the way towards a many-body picture '®^ of structural 
instabilities in low-dimensional electronic as well as magnetic chains. The key asset of many- 
body approaches resides in the possibility to map the ID electron-phonon problem onto a 
ID interacting electron gas^"^^^ with retarded electron-electron interactions, so that well- 
developed quantum-field theoretical concepts for low-dimensional electronic liquids can also 
apply in these specific contexts.'''' ^'^^^ In particular, this mapping revealed the inconsistency 
of the quasiparticle picture inherent to one-electron perturbative approaches, and permitted 
a better understanding of several experimental observations one must fascinating being the 
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non-universal character of critical exponents of correlation functions. ^"^^' ^^' 
The present study aims at investigating the charge-density-wave(CDW) instability for the 
molecular-crystal model with quartic anharmonic phonons, with emphasis on the contribu- 
tion of phonon anharmonicity to the transition temperature. A key step in this respect will 
be an exact formulation of the effective two-electron action for this model which is valid at 
any arbitrary regime(i.e. weak or strong) of crystal anharmonicity. Instructively, the issue of 
the interplay of electrons and anharmonic phonons has been considered recently ^"^' ^^"^"^^ in 
attempts to explain unsual local structure effects on thermal and spectral properties of cer- 
tain materials.^' ^^'"'^^' In particular, refs.^^'"*^^' have addressed the specific problem of phonon 
anharmonicity and its effects on the superconducting and CDW transition temperatures. 
Thus, following a self-consistent Gaussian decoupling of the quartic anharmonic phonon, it 
has been established that the electron-anharmonic phonon model could be reduced to an 
effective electron-harmonic phonon one without the particle-hole symmetry, but with a per- 
sistent finite transition temperature at half filling of the electronic band. On the basis of this 
self-consistent Gaussian treatment, the authors found that the phonon anharmonicity had 
no qualitative effect on electronic instabilities in anharmonic crystals except breaking the 
particle-hole symmetry. 

Here we wish to extend results of refs.^^' ^^^ for the CDW instability, to materials with relatively 
strong quartic anharmonicities. In this purpose, we shall establish that by proper account of 
phonons via non-Gaussian partial functional integrations^^' ^^' of phonon degrees of freedom, 
one can obtain the exact retarded two-electron action for the anharmonic molecular-crystal 
Hamiltonian. By use of appropriate Hubbard-Stratonovich fields, an effective field theory for 
the CDW instability will be constructed and the transition temperature derived. 

2. The anharmonic-phonon-induced two-electron action 

The standard Holstein model^' describes a ID electron gas coupled to a molecular lattice 
with quasi-optical harmonic vibrations. For crystal lattices dominated by quartic anharmonic 
vibrations, this model can be extended leading to an anharmonic Holstein model which, at 
half-filling of the electronic band, reads: 



H 



^ 'Pf Mu^'o 2 . 4 

— — H r + A r ■ 



E 



N N 

— t 2_^ '^j,a- Cji^a + O! 2_^ ^j Cj^a Cj,a, (Ij 

where rj is the displacement of the j molecule of mass M and characteristic frequency cOq 
relative to its equilibrium, Pj is the momentum conjugate to rj and A is the coefficient of 
anharmonicity. Cj ^{cj^a) creates(annihilates) an electron state of spin a at site j and t is the 
uniform transfer integral between nearest-neighbour sites j and j' = j it 1. Lastly, a is the 
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electron-phonon interaction coefficient. 

Tlie Holstein model fias a long and rich history in the literature of polarons in ID and 2D 
materials, high-temperature cuprate superconductors and colossal magnetoresistance man- 
ganites^^"^'^' being two classes of materials that attracted much recent attention because of a 
sizable signature of polaronic effects in their heat and charge transport parameters. As for elec- 
tronic instabilities, this model has been shown to possess a rich phase diagram covering CDW 
and superconductivity as observed in Bai-xKxBiO^, a material that exhibits a CDW-driven 
metal-insulator transition and which upon doping becomes superconducting. From the stand- 
point of fundamental physics, the Holstein model is particularly interesting as it represents 
one among the few rare models in which a net departure from Fermi liquid behaviour occurs. 
Indeed, there the electronic spectral density looses the quasiparticle peak and is swamped 
by an incoherent background with increasing electron-phonon coupling. ^^'^^' In general these 
properties can be unveiled either analytically as in the case of the ID Holstein model with 
a linearized electronic dispersion which is exactly solvable, or via quantum Monte-Carlo and 
diagrammatic approaches based on Migdal-Eliashberg theory as well established^^"^^' for the 
2D Holstein model. At half filling of the electronic band and strong electron-phonon cou- 
pling, the model gives rise to an effective Hubbard model so that CDW and on-site s-wave 
superconductivity coexist. Away from half filling or with decreasing electron-phonon coupling 
but high phonon frequencies, the superconducting phase dominates. In recent years there has 
been much interest to the Holstein-Hubbard model^^"^''^ for which the competition between 
phonon-mediated and direct electronic correlations, permits the account of new groundstates 
in the phase diagram of the bare Holstein model including spin-density-wave and d-wave su- 
perconducting phases. Remarkably, the functional formalism involving partial integrations of 
phonons has demonstrated great effectiveness in this last context by furnishing the retarded 
attractive two-electron Hamiltonian competing with Hubbard's Coulomb repulsion.^^'^^'^^' 
In this work we address the issue of a possible exact formulation of the phonon-mediated 
two-electron action for the Holstein model with quartic anharmonic phonons (1), given its 
importance that extends far beyond the specific context of CDW instability considered in the 
present context if looked out from the more global standpoint of the effective-field theory for 
electronic instabilities within the framework of the Holstein-Hubbard model. ^"^^^ Partially in- 
tegrating out lattice degrees of freedom in the total Hamiltonian (1), we arrive at the following 
partition function in the Matsubara space and field-theory representation:^^) 



where 



zO=/n^r(g,r)e^?« (3) 



9,T 
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is the bare lattice partition function corresponding to the anharmonic action: 



V / dTr{qi,T)r{q2,T)r{qs,T)r{q4,T) 



gi,<?2,<?3,<?4 

(4) 

L = Na being the total length of the system. In formula (4), r is the bosonic degree of freedom 
in the field configuration for the classical displacement r i.e. 

r{x,T) = —=y2r{q,T)exp{iqx), x = j a, (5) 

while 'Do(r), the free-phonon propagator, is derived from the general expression of the retarded 
one-phonon propagator i.e. 

Vx{T) = -<Trr{q,T)r{-q,0)>x, (6) 

in the Gaussian limit A — t- 0. In the Matsubara space T>o{t) becomes: 

Vo{uJrn) = 2~~^ 3' C^) 



with ujm = 2TTmf3 ^ the Matsubara frequency. Lastly, 

'5° = E / dTrkA^)G,\k,T)^PkAr) 



(8) 



is the free-electron action where 



Go{k,T) = -d/dT-VFk (9) 

is the free-electron propagator for an electronic dispersion linearized around the Fermi wavec- 
tor kp{vp being the Fermi velocity), and ip is the usual Grassmann field. 
The quantity in (2) of primary interest to us is the phonon-mediated effective two-electron 
action, i.e.: 

Si{i^,i^*) = In [< e5e-KV',r;r) > 1 ^ (10) 

in which 

S,^liij,r; r) = ^Y.J dTrk+gA^)i^kAr)riq,r) (11) 

k,q,(T 

is the electron-lattice interaction action with averages over the anharmonic lattice action. 
Since the lattice action is anharmonic in the phonon fields, a Gaussian treatment of averages 
over phonons requires that the anharmonic coefficient is very weak compared to the harmonic 
stifness Mujq.^'^' However, if the two coupling coefficients are strongs the above perturbation 
consideration will inevitably result in a crude approximation of the problem. Actually, only 
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by exact integration of phonon degrees of freedom can one gain full understanding of both 
qualitative and quantitative contributions of the phonon anharmonicity to the retardation 
process in the resulting effective two-electron action. In this last purpose, assume for now 
that the quadratic mean over boson fields in (lO)-(ll) can be calculated exactly and leads to 
an effective interaction action which we formally represent as: 

Si = -^ E / d^' I '^^axir - t') rk+,A^)rk'-q,A^')^k'Ar')i>kAr)- (12) 

k,k',q,a,a' "^^ -^ ^ 

In the last formula the phonon-mediated electron-electron interaction matrix has been defined 
as: 

gx{T-T') = a^Vx{T-T'), (13) 

where the quantity T)\{t — t') is the exact one-phonon propagator in the presence of quartic 
anharmonicity. Following conventional notations, we define the Matsubara Fourier transform 
of V\{t — t') as: 

oo 

Vx{T-r')=r' E ^xAm)e-'^-^^-^'\ (14) 



m=— oo 



and with the help of a standard non-Gaussian integral, ^^'^^' the anharmonic one-phonon 
propagator Vx{uJm) is exactly found: 

'"^-3/2 [ 



1 ^-3/2[;^^0"'(^™)] 
^A(Wm) = 



y2A ^-1/2 [^^0 H^ni)] ' 

A = XT/L, (15) 

with Dp{y) the Parabolic Cylinder function of order p and argument y.'^^' With help of the 
explicit formula of T)\[ujra) just derived, we transform the two-electron action Sj in the Mat- 
subara Fourier space obtaining: 

T 

k,k' ,q,cr,cr' 

where k = {k,u}n), q = {q,u}m) and the phonon-mediated electron-electron interaction matrix 
now expressed in Matsubara Fourier space reads: 

1 ^-3/2[;^^0"'(^™)] 

axArn, T) = -A —= _i ^^ -7. (17) 

V2A ^-l/2[^=f^0 i^^rn)] 
Y zA 

To see what formula (17) becomes in the harmonic-phonon regime, remark that when A — )■ 
the argument of the Parabolic Cylinder function becomes infinitely large. Thus, we can readily 
use the following asymptotic expansion of the Parabolic Cylinder function:'^^' 

D-n-l/2{y) CC y-"-V2[i + o(lMexp(-^), 



5/12 



J. Phys. Soc. Jpn, 
2.0 



Full Paper 




Fig. 1. The reduced one-phonon propagator versus iwm/wo, for Ae// = 0, 0.2, 0.5, 1, 2, 5. The diverg- 
ing curve corresponding to Ae// = 0, coincides exactly with the reduced free-phonon propagator 
I?o(wm)/2?o(0). 



oo, 



to rewrite the anharmonic one-phonon propagator (15) as: 



(18) 



(19) 



The weU-known harmonic result for the phonon-mediated electron-electron interaction matrix 
then follows: 

= -a'/Miul + Lol). (20) 

In the classical regime i.e. m = 0, go{0) is negative implying a retarded attractive electron- 
electron interaction. 

In the general case when the anharmonicity coefficient A is of arbitrary strength, the physics 
behind formula (17) is not easy to capture from simple analytical considerations. Proceeding 
by numerical analysis, it is useful for a better understanding of the role played by anhar- 
monic phonons in the two-electron correlations reflected by the interaction action (16)-(17), 
to first examine the salient features of the one-phonon propagator (15). On figure 1, we plot 
Dx{ujm)/DQ{0) versus the reduced real frequency iujm/i^o for different values of the effective 
anharmonicity coefficient Ag// = X/{Mio'^)'^. When Ag// = 0, the one-phonon propagator is 
singular at iuJrn = ^o in agreement with established behaviour of the free-phonon propaga- 
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Fig. 2. Variation of the reduced effective electron-electron interaction gx{Q,T)/ go{i)) as a function 
of the effective anharmonicity coefficient Ag//, for iLOm = 0, 0.25 wq, 0.5 wq, 0.75 wq and a;o(froin 
bottom to top curves). 

tor. However, when the effective anharmonicity coefficient becomes nonzero the singularity is 
completely suppressed leading to optical phonon states of finite lifetimes. ^^^ As we increase 
Ae//, the finite tail of double peaks in the phonon propagator at the characteristic frequency 
Wo decreases and vanishes at relatively strong anharmonicities relative to the characteristic 
harmonic-phonon stifness Mwq. Figure (2) displays the resulting curves showing the variation 
of gx{ujm,T)/gQ{0) as a function of Ae//, for iujn = 0, 0.25ujq, O.Swq, O.TSwq and a;o(from 
bottom to the top curves). 

Instructively, in addition to the pure harmonic regime consistently described for A = 0, for- 
mula (17) also suggests two other distinct physical contexts involving a vanishing effective 
electron-electron interaction gx{ujm, T). They are L — )• oo but T finite, and T = but L finite. 
These two characteristic features, which are proper to the exact treatment of the anharmonic 
phonons, give evidence of the contribution of material sizes to electronic correlations mediated 
by anharmonic lattice vibrations. 

3. Anharmonic-phonon-induced CDW instability 

In this section we shall consider the relevant problem of CDW instability which has been 
widely studied within the framework of the Holstein model with harmonic lattice vibrations, 
and examine the effect of phonon anharmonicity on the transition temperature. Proceeding 
with, we follow the effective-field approach which involves constructing an effective Landau- 
Ginzburg(LG) functional associate with CDW fluctuations governed by the effective two- 
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electron action (16). First, we define the CDW operator: 

0{q, r) = -^ V r{.K r)Hk + q, r) (21) 

which helps rewrite the effective interaction action in a suitable form i.e.: 

Si = -\y^ I dr' f dTgx{T-T')Oiq,T)0{-q,T'). (22) 

Next, introducing appropriate Hubbard-Stratonovich fields conjugate to the above CDW opra- 
tors, the total partition function (2) becomes: 

X exp{-Y, dr' dT4>{q,T)[2g^\T - T'm-q,T')} 
q Jo Jo 

X exp{-T dT[4>{q,T)0{-q,T)+4>{-q,r)0{q,T)]}. (23) 

Upon total elimination of electronic degrees of freedom and use of the linked cluster theorem, 
the last expression reduces to: 

Z = zOe-^-^[^l, 

Z' = Z^Zl (24) 

where 

/3-^M = Y. [^9l\i^m,T) - x{i^m,T)] |0(cz,u;™)|2 + ... (25) 

is the LG functional in which the CDW-CDW correlation function: 

O'T' 

x{oJm.,T) = — — ^Go(/c,a;„)Go(A;- 2fcF,a;„ -cjm), Wn = (2n + l)/3"\ (26) 

k.UJn 



has been defined and is explicitely obtained for the nearly- free electron system as^ ' ^^ 



x{uJm,T) = Nf 



ln^^ + ^m-^f4+"^^ 



(27) 



T " \2J \2 4ttT , 

with Np = l/2TrvF, ip the Digamma function and Ep = vpkp the Fermi energy. In the 
classical regime(i.e. oJm = 0), the soft-mode condition requires the vanishing of coefficient of 
the Gaussian term in the LG functional (25). To extract the critical temperature from this 
condition it is convenient to introduce T^ for the harmonic-phonon Holstein model i.e.: 

TO = 1.13^p exp 4, 90 = NpgoiO). (28) 

With the help of this parameter we can define a reduced CDW transition temperature 
tc = T^/Tj}, and by numerical root finding we extract the Ag// — tc phase diagram shown 
in figure 3. More precisely, on the figure t^ is plotted versus Ae// for four representative values 
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Aeff 



Fig. 3. The reduced CDW critical temperarture versus the effective anharmonicity coefhcient Xeff, 
for a'^/MujQ = 0.05, 0.1, 0.5 and l(from bottom to top curves). 



of the bare electron-phonon interaction coefficient go{0), assumed to reflect distinct regimes 
of the electron-phonon interaction. A most remarkable fact emerging from curves is that tc 
is alvi^ays a decreasing function of Ae//, irrespective of the strength of the electron-phonon 
interaction. Equally remarkable is the reduction of the CDW transition temperature which 
is manifestly as strong as the electron-phonon interaction is weak. Quite strickingly no re- 
entrance phase appears for all the representative values of the electron-phonon interaction 
considered, contrasting with predictions of the self-consistent Gaussian phonon approxima- 

4. conclusions 

We obtained the expression of the phonon-mediated two-electron action for the molecular- 
crystal model with quartic anharmonic phonons, valid at any order of anharmonicity strength. 
To achieve our goal, we followed a non-perturbative approach based on non-Gaussian partial 
functional integrations of anharmonic phonon degrees of freedom. As a practical illustration 
of the advantage for integrating out exactly the anharmonic phonons we considered the CDW 
problem, and shown that the account of quartic anharmonic phonons in the Holstein model 
gives rise to a picture of the transition dominated by the interplay of material's local structure 
including the volume. However, for fixed volume it emerged that the suppression of the CDW 
transition was as strong as the electron-phonon interaction was weak. No re-entrant process 
was observed in the weak, medium and strong electron-phonon coupling regimes in contradic- 
tion with a recent prediction based on the self-consistent Gaussian approximation.^^' ^-^ 
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From a general standpoint, anharmonic lattices display a wealth of rich properties that are 
distinct from their harmonic counterparts as established at length in past theoretical as well 
as experimental researches. Of fundamental interest, multi-phonon correlations lead to un- 
usual phenomena such as biphonon boundstates, a broadening of one-phonon modes which is 
not a simple shift of mode frequencies in Raman spectra and the appearance of a distinctive 
structure of the phonon spectrum dominated by modes of finite lifetimes. Unveiling these 
properties requires a non perturbative treatment of anharmonic lattice vibrations in view of 
full understanding of their contributions to structural instabilities of anharmonic crystals. 
On the other hand, the importance of anharmonicity in electronic phase instabilities of low- 
dimensional materials has been demonstated in various experimental probes on lattice and 
electronic structures, including extended x-ray-absorption fine structures, neutron scatter- 
ing and optical techniques. The lattice anharmonicity is thus observed as being particularly 
pronounced in planar structures such as cuprate oxydes and colossal manganites where it 
governs both intra and inter planar charge transfers. In most of these materials, the plane 
chain structure provides a suitable scenario for charge-transfer excitations owing to a possible 
reduction of the Coulomb repulsion due namely to interplane interactions, by anharmonic 
vibrations of ions along the chain(as in cuprate oxydes, see e.g.^'). In this context the phonon 
anharmonicity can be so strong that the shift of the one-phonon frequency in the infrared Ra- 
man spectrum predicted within the self-consistent Gaussian renormalization approach, is not 
meaningful for true characterization of the temperature dependence of the effective charge- 
transfer interaction, the contribution of the local structure to structural instabilities and so 
on^'^^-* as observed experimentally in real anharmonic materials. 

Note to end that in recent years, there has been a great deal of interest^"^"'^'"'^'^-' to the interplay 
of phonon-mediated attraction and Coulomb repulsion between electrons in Luttinger liquids. 
In particular, the Holstein-Hubbard model has been shown to display a phase diagram char- 
acterized by a great variety of groundstates including Mott insulator, charge-density-wave, 
spin-density-wave, s-wave and d-wave superconducting groundstates. For anharmonic mate- 
rials belonging to this class, an exact formulation of the anharmonic phonon contribution to 
electronic correlations represents a relevant prerequisite towards a best account of specific 
features of their phase diagrams. 
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